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We investigate waveguiding of ultraslow light pulses in an atomic Bose- 
Einstein condensate. We show that under the conditions of off-resonant 
electromagnetically induced transparency, waveguiding with a few ultraslow 
modes can be realized. The number of modes that can be supported by 
the condensate can be controlled by means of experimentally accessible 
parameters. Propagation constants and the mode conditions are determined 
analytically using a WKB analysis. Mode profiles are found numerically. © 
2008 Optical Society of America 
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An optical pulse can propagate at ultraslow speeds through an atomic Bose-Einstein con- 
densate (BEC) 1 under electromagnetically induced transparency (EIT) 2 conditions. Such 
pulses can be utilized for storage of coherent optical information. 3 ' 4 Most of the early the- 
oretical investigations assume one-dimensional light propagation, 5 except for a few notable 
studies. 6-8 Confinement of optical pulses in all three directions has been examined quite 
recently. 9 Transverse intensity profile of a focused control pulse translates to an effective, 
spatially varying refractive index for the probe pulse, which can be used for waveguiding 
and confinement purposes. 9 Controlled spatial variation of the refractive index profile was 
first suggested and used for creation of quantum shock waves with ultra-compressed slow 
light pulses. 6 The shock wave induced snake instability, depending on the transverse density 
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variations of the BEC, was observed. 6 More recently, transverse effects on slow light have 
been studied numerically. 10 It has been found that, using waveguiding by control pulses with 
Gaussian shaped transverse profiles, smaller group speeds can be achieved for higher order 
modes of the probe pulse. 10 

Possibility of slow light waveguiding in ultracold atomic medium has been demonstrated 
in a recent experiment. 11 Decrease of refractive index away from the optical axis and the 
corresponding graded index optical waveguiding are consequences of Gaussian density profile 
of the ultracold atomic cloud, strongly confined in transverse dimensions. 9 The experiment 
uses recoil induced resonance (RIR) method to achieve slow light (c ~ 1500 m/s), due to the 
strong dispersion in the high gain regime. 9 The ultracold atomic medium with RIR scheme 
operates as a graded index optical waveguide in the strongly guided regime due to the large 
core radius (~ 200 /mi) and high refractive index contrast. 11 

In this letter, we investigate ultraslow waveguiding by BEC under EIT conditions. By con- 
sidering slightly off-resonant EIT pulses propagating in a BEC, tightly trapped in transverse 
dimensions, possibility of weakly guided regime with few modes is shown. Single mode con- 
dition is established. These regimes are promising for guided nonlinear optical phenomena 
in ultracold matter. 11 Multiple modes at low temperatures are found. These results might 
be useful to design spatially controllable higher capacity optical memories. 10 

At low temperatures a Bose gas can be considered as a condensed cloud in a thermal 
gas background. Atomic number density profile of such a system can be described by 12 
p(r) = p c (r) + Pth(^), where p c {r) = [(p — V(r))/U ]Q(p — V(r)) is the number density 
of the condensed atoms and p t h is the density of the noncondensed ideal Bose gas. Here 
Uq = 4irh 2 a s /m; m is atomic mass; a s is the atomic s-wave scattering length. 0(.) is the 
Heaviside step function and Tq is the critical temperature. The external trapping potential 
is V(r) = (m/2)(u; 2 r 2 + oj 2 z 2 ) with u r ,uj z are trap frequencies for the radial and axial 
directions, respectively. At temperatures below T c , p, the chemical potential is determined by 
/i(T) = p TF (N /N) 2 ^, where Ptf is the chemical potential evaluated under Thomas-Fermi 



approximation, ptf = ((^t) /2)(15iVa s / ah) 2 ^ , with cu t = (lUzLU 2 ) 1 / 3 and ah = yh/(cu z ujf.) 1 / 3 , 

the average harmonic oscillator length scale. The condensate fraction is given by Nq/N = 
1 — x 3 — s(C(2)/^(3))x 2 (l — x 3 ) 2 / 5 , with x = T/T c , and £ is the Riemann-Zeta function. The 
scaling parameter s is given by s = ptfI^bTc — (l/2)C(3) 1 ^ 3 (15iV 1 / 6 a s /a/ l ) 2/ ' 5 . 
EIT susceptibility 2 for BEC of atomic density p can be expressed as \ = pxi with 



where A = uo — ujq is the detuning of the probe field frequency uj from the atomic resonance 
uj . Q c is the Rabi frequency of the control field; p is the dipole matrix element for the probe 
transition. T 2 and T 3 denote the dephasing rates of the atomic coherence. We consider a 




Xi = 



\pI i(-iA + r 2 /2) 

e h [(r 2 /2 - iA)(r 3 /2 - iA) + fi 2 /4] ' 
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BEC of 23 Na atoms with parameters, 1 N = 8.3 x 10 6 , u r = 2ix x 69 Hz, uj z = 2ir x 21 Hz, 
r 3 = 0.57, 7 = 2vr x 10.01 MHz, and T 2 = 2vr x 10 3 Hz. We take fi c = 2.57 and A = -O.I7 
so that x' = 0-04 and x" — 0.0006 at T = 42 nK. Here, x' and x" are the real and imaginary 



parts of X) respectively. Neglecting x" ■> the refractive index becomes n = a/1 + x' ■ in the 
thermal cloud, x' ^ 1 so that an approximate refractive index profile in the radial direction 
can be written as 

, . / m[l - Am 2 ] 1 / 2 r , . 

n ( r ) = s ^ „ , (2) 

II r > R 



where A = 1 — 1/n, 2 and i? = y2/i(T)/muj^. In general ni is the index along the center 
line (z axis) of the condensate. n\ is different than 1 only in the vicinity of the center of 
the condensate, where it varies slowly over the order of the optical wavelength. It can be 
determined by n\ = (l+fiXi/Uo) 1 ' 2 - In the radial direction, BEC density profile is translated 
to a refractive index analogous to a graded index fiber. Thermal gas background plays the 
role of the fiber coating while the condensate acts as the core. 

We use the cylindrical coordinates as the refractive index n(r) is axially symmetric. Small 
index difference between the condensed and noncondensed clouds results in weak guiding 
of probe pulse, for which the normal modes are the linearly polarized modes (LP modes) 13 
which are determined by the Helmholtz radial equation 

d 2 1 d 



^^-^ + ^ 2 (r)-/3 2 -/ 2 /r 2 
ar z r ar 



^(r) = 0. (3) 



Here ko = 00 /c, I = 0, 1, 2, (3, and ip are defined for the transverse field of the LP modes 
E t = ip(r) exp [i(/0 + ut — f3z))}. Employing Wentzel-Kramers-Brillouin (WKB) method, the 
characteristic equation is found to be 

ri 7 2 1 

[kln\r) - A 2 m - ^] 1/2 dr = (m + -)vr, (4) 
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where 77 2 are the zeros of the integrand and m = 0, 1, 2, 3, .... Solving Eq. [4] for fl, we find 

A m = mfco[i-^^] 1/2 , (5) 

riik Q R 

where K = (/+2m+l). The propagation constant /3i m , depends on the frequency of the probe 
pulse. This frequency dependence determines the dispersion of the waveguide. Propagation 
constant can also be defined in terms of an effective mode index ri[ m such that (3i m = rii m k = 
ni m {n u u)u/c. 

The group velocity v g is defined as l/v g = d/?/da>, which leads to 

1 d(3 uj ( driim dni drii m \ n lm 



v q du c \ drii du du ) c 



(6) 
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Here, the terms in the parenthesis are due to the material and waveguide dispersion, respec- 
tively. Under EIT conditions, material dispersion dominates. Furhtermore, n\ m ~ 1 so that 
we can write 

The group velocity at T = 42 nK as a function of K is given in Fig{TJ We calculate v g = 30.8 
m/s, for the LP o mode. Group velocity increases with the mode number K. The higher 
order modes are away from the high density region of the condensed cloud and move faster. 
This behavior reduces modal dispersion and can be beneficial to preserve shape of the probe 
pulse for coherent optical information applications in BECs. 

The dimensionless normalized frequency V is defined as V = nik R\^A = k R(n1 — 
l) 1 / 2 . It determines the number of modes that can be supported by the atomic cloud. We 
find a confined mode description via k$ < (3 < riik so that K < V/2. A condensate 
supports only single LP o mode for 2 < V < 4. For 2 < V < 6, it could support just 
two modes, LP o and LP 10 . V can be tuned by various experimental parameters as V = 
ko[i(2x'i/muJr u o) 1 ^ 2 - Multiple-mode support can be perhaps most conveniently tuned by T. 
Temperature dependence of V is plotted in FigfSJ For the parameters used in Fig. [TJ at T = 42 
nK it starts from V = 45 and decreases down to a lowest value of V ~ 7.3 at T = 398.7 nK. 
As temperature increases, the condensed cloud shrinks and the index contrast diminish so 
that the BEC can no longer support multiple mode propagation. 

More accurate results may be obtained using a numerical scheme. For that purpose it is 
assumed that the index of refraction increases incrementally from the edge to the center of 
the core. If the number of increments is high enough, this model will be a faithful represen- 
tation of the conditions prevailing in the inhomogeneous core. Within each increment the 
refractive index is constant and we can use exact solutions of the wave equation in terms of 
Bessel functions. Electromagnetic boundary conditions provide a recurrence relation between 
successive solutions as well as an equation for the propagation constant. We have compared 
the results of our numerical calculation with these of WKB calculations and found excel- 
lent agreement. Typical results of our numerical simulations are given in Figj3l where the 
modes LP o and LP 10 are shown. The modes are localized in the condensate core of radius 
R — 21.1 /im. Beyond R, in the thermal component, the modes are evanescent. Such modes 
can be addressed by specifically constructed transverse profiles of the probe pulses. This can 
be exploited to control of capacity in the transverse direction for storage of coherent optical 
information in different mode patterns. These results can be further combined with the lon- 
gitudinal control of optical information storage capacity where a particular EIT scheme, in 
which the control field is perpendicular to the probe field, is employed. 6 

Summarizing, we have examined optical waveguiding of ultraslow pulses in a Bose-Einstein 
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condensate under EIT conditions. Off-resonant EIT scheme has been considered to permit 
multi-mode ultraslow light propagation. Propagation constant is calculated analytically us- 
ing WKB approximation. Mode profiles are determined numerically. Single and two mode 
conditions are established in terms of experimentally accessible parameters. Temperature de- 
pendence of the number of modes that can be supported by BEC is presented. Such ultraslow 
modes may be useful for realizing spatially controllable storage of coherent optical informa- 
tion with higher capacity in BECs, in the transverse direction, complementing longitudinal 
control of information storage capacity. 6 

O.E.M. acknowledges support by TUBA-GEBIP Award. D.T. is supported by Istanbul 
Technical University Foundation (ITU BAP). 
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List of Figures 

1 Dependence of the group velocity on K = (/ + 2m + 1) at T = 42 nK (T c = 424 
nK). Parameters are given in the text. For the single mode condition group 
velocity is 30.8 m/s 

2 Temperature dependence of V is plotted from T = 42 nK to T = 398.7 nK 

(T c = 424 nK). V is related to K via K < V/2 | 

3 Radial profiles of transverse modes LP o (solid line) and LPi (dotted line) 
confined in a BEC of 23 Na atoms with parameters given in the text 
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Fig. 1. Dependence of the group velocity on K — (I + 2m + 1) at T = 42 nK (T c = 424 nK). 
Parameters are given in the text. For the single mode condition group velocity is 30.8 m/s. 



7 



50 




40 


— 




30 






20 






10 














0.2 0.4 0.6 0.8 1 



T/T c 

Fig. 2. Temperature dependence of V is plotted from T = 42 nK to T = 398.7 nK (T c = 424 
nK). V is related to K via K < V/2. 
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Fig. 3. Radial profiles of transverse modes LPoo (solid line) and LPio (dotted line) confined 
in a BEC of 23 Na atoms with parameters given in the text. 
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